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We study the emergent proximity effect in mesoscopic circuits that involve a conventional su- 
perconductor and an unconventional pnictide superconductor separated by a diffusive normal or 
ferromagnetic wire. The focus is placed on revealing signatures of the proposed s" 1 state of pnic- 
tides from the proximity-induced density of states and Josephson current. We find analytically a 
universal result for the density of states that exhibits both the Thouless gap at low energies, and 
peculiar features near the superconducting gap edges at higher energies. The latter may be used to 
discriminate between s H and s ++ symmetry scenarios in scanning tunneling spectroscopy experi- 
ments. We also calculate Josephson current-phase relationships for different junction configurations, 
which are found to display robust 0-tt transitions for a wide range of parameters. 
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I. INTRODUCTION 

The origin of unconventional superconductivity in fer- 
ropnictide compounds, their phase diagram and symme- 
try of the underlying order parameter are topics that 
attract considerable interest in recent years, see Refs. Q3- 
[3] for reviews. Superconductivity in pnictides emerges in 
close proximity to an antifcrromagnetically ordered state, 
and the critical te mpe rature T c has a dome-shaped de- 
pendence on doping^ similar to that in cuprates. Due to 
their multiband electronic structure with multiple Fermi 
surfaces and delicate interplay of interactions in different 
channels a number of possibilities for electron ordering 
are possibleP^ Structural transition, competing or co- 
existing magnetic spin-density-wave (SDW) and super- 
conducting (SC) orders are being examples. The latter 
may be in the form of the conventional s ++ -wave state 
that has s-wave symmetry in the Brillouin zone and gaps 
of the same sign on electron and holes Fermi surfaces. 
Alternatively, SC order may appear in the form of an 
extended s~* state that looks as s-wave from a symme- 
try point of view but has opposite signs of the gaps on 
different sheets of the Fermi surfacesPtETJ There may be 
scenarios of several SC states with the nodes in the SC 
gap, of both s-wave and d-wave symmetries JA2HH1 

Experimentally, the most convincing support in favor 
of unconventional symmetry of pnictides is given so far 
by the observed spin resonance below T c in inelastic neu- 
tron scattering measurements on K-doped BaFeAsP^ In 
all the materials studied, the resonance occurs at the an- 
tifcrromagnetic wave vector Q of the parent compound. 
It is thought to be a triplet excitation of the singlet 
Cooper pairs, implying a superconducting order param- 
eter that satisfies A(k + Q) = — A(fc), which indicates 
either s H or d-wave cases. Fabricated c-axis Josephson 
junctions of this material and ordinary superconductor 
are suggestive of an s-wave state, but not providing un- 
ambiguous evidence for the s H state itself ™ In addition, 
a&-corner-j unction experiments with Co-doped BaFeAs 
seems to eliminate the option of d-wave pairing.^ Other 
notable experiments providing substantial but still indi- 



rect support of s H state include quasiparticle interfer- 
ence in magnetic field probed by scanning tunneling mi- 
croscopy^ and observation of half-integer flux-quantum 
jumps through the loop formed by niobium and polycrys- 
talline iron-pnictide sample Finally, there is a growing 
number of low-temperatures studies addressing thermo- 
dynamics and transport properties of pnictides, however 
it is usually hard to deduce underlying symmetry of a 
superconductor from such data. For fully gapped s-wave 
state, one expects to see exponentially suppressed quasi- 
particle response and power-law in temperature for the 
d-wave state with the nodes. The possible ambiguity in 
interpretation of data stems from the fact that accidental 
nodes on the Fermi surface or impurity-induced subgap 
states may easily alter low-temperature behavior of, for 
example, heat capacity or London penetration depth. 

It is widely agreed that a decisive experiment should 
involve a phase sensitive probe such a s the Joseph- 
son effect. Although original proposal d 20 * 21 * followed 
mostly immediately after the s" 1 candidate symme- 
try was introduced, no such direct measurements of the 
current-phase relationship have been performed so far for 
pnictide-based Josephson junctions. Nevertheless, this 
inspired a lot of theoretical efforts in finding simpler ge- 
ometries or alternative signatures of pairing state in 
proximity circuits with pnictides and conventional super- 
conductors! 2 ^"^ A particularly interesting recent conclu- 
sior l 24 1 29 1 is that tunneling spectra of weakly coupled s- 
s H bi-layers exhibit distinct features characteristic only 
to sign-changing symmetry of the gap. Physically, the 
effect comes from the frustration in the junction since 
the gap of an ordinary superconductor tends to align 
with one of the gaps of s H superconductor, and thus 
becomes in the conflict with the other band experiencing 
the anti-proximity effect. 



II. THEORETICAL FRAMEWORK 

In this work, we consider plethora of effects in 
superconductor-metal-superconductor junctions where 
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one or both superconductors are assumed to have s H 
symmetry. The metal is either normal or ferromagnetic 
diffusive wire. In the context of the Josephson effect, 
in such structures, we find various current-phase rela- 
tionships whose shapes depend on the relation between 
the wire length and superconducting coherence length, 
and boundary transparency. The generic feature is non- 
analytical behavior of the current near phase tt, which 
corresponds to the closing of the proximity-induced gap 
in the wire, and robust 0-tt oscillations even without fer- 
romagnets. In the context of the proximity-induced den- 
sity of states (DOS) in the wire, we identify fingerprints 
of the s" 1 symmetry, which is thus not only unique to 
frustrated s-s H bi-layers. 

We build our calculations based on the Usadel equa- 
tions^! and accompanying Kupriyanov-Lukichev bound- 
ary conditions^ modified for a multiband caseP^ This 
quasiclassical theory captures all the essential features 
and full complexity of the proximity effect. Adopting 
angular parametrization for the normal and anomalous 
quasiclassical Green's functions^ as G = cos 9(u>,x) and 
F = sin6*(oj, x)e %x ^ w,x \ Usadel equations take the form 



dl6 -(2u/e Th ) sin 6 = {d xX f 



sin o cos ( 



^(sin 2 ^ X ) =0, 



(la) 
(lb) 

where uj = (2n + 1)itT is Matsubara frequency, Etk — 
D/L 2 is the Thouless energy for the wire of length L, 
and D is the diffusion coefficient. Spatial derivatives are 
taken with respect to the dimensionless coordinate x — > 
x/L and we assume quasi-one-dimensional geometry. At 
the interface, we have two boundary conditions: 

Ju> = 2 (<5A/7A)sm0Bsm0 sA sinV>, (2a) 



A=l,2 



s-^ cos 9 sX sin 9 B 
2 > [b\ tan 6 S \ cot 9 b cos ip — 1J 

A=l,2 



7a 



(2b) 

Here, J w denotes the first integral of Eq. (lb), index A 
labels different bands, ip = (4> — 2xb)/2 and <t> stands for 
the global superconducting phase difference across the 
junction, while factors 8\ = ±1 account for the rela- 
tive shifts of phases between the bands, and finally pa- 
rameters 7a represent dimensionless interface resistances. 
We also used notations 9b — 9 (lu, ±1/2) and similar 
for XBi and introduced Green's functions of a super- 

uj 2 and 



conductor in the bulk: sin6> s A = |Aa|/\/|A;, 



cos 9* 



uj 2 , with A\ being corresponding 
gaps. Having solved Usadel equations, one can find a 
density of states 



N(e,x)/N Q = Re[cos(w, x)] u 



(3) 



upon analytical continuation to real energies, and a 
Josephson current-phase relationship 



eI{4>)RN — 2-7rT^ Ju = J tanh -^liaJ e de 



(4) 



upon summation over Matsubara frequencies, where Nq 
is bare density of states in a metal and Rn = L/e 2 DN S 
is normal state wire resistance of cross-section area S. 



III. DENSITY OF STATES 

Consider a symmetric s\n\s junction. In the absence 
of superconducting phase difference between the leads we 
have x = and the whole system of equations simplifies 
to one: 



dl6- (2oj/e Th ) sin = 0, 



(5) 



which has to be solved for x G [—1/2,1/2]. Since this 
is the same equation as for the nonlinear pendulum it 
can be integrated exactly in terms of the Jacobi elliptic 
functions. Indeed, the above equation has a simple first 
integral: 



(8 X 9) 2 = (4u/ez7»)[caB0 o -cos0], 



(G) 



where integration constant 9q = 9(lu,0) was chosen to 
be at the middle of the wire due to obvious symmetry 
reasons. To perform a subsequent second integration, we 
change variables as 



cos f = 



2m cos 2 </ 
1 — m sin 



and find 



2w 



1, m = cos 2 (0 o /2), (7) 



(8) 



m sin 



which is a tabulated integral. Finally, using the Jacobi 
elliptic functions cn and dn, the solution appears in the 
form 



cos[9(lo,x)/2] = cos(6» /2) 



cn^u, m) 
dn(u, to) ' 



2uj 

£Th 



It is important to keep in mind that the modulus of the 
Jacobi functions is actually an energy-dependent function 
to(oj), see Eq. ( 7). By using now Eq. ^ in the boundary 



condition Eq. 



, one finds a closed algebraic equation 



for the unknown integration coefficient in the form 



Ub 



sn(u^,TO) cn(us,TO) 
dn(uB, m) 



I- 



\J to(1 — to) 
where ub = \/w/2sTin and 

G(w) = cos s xhx = ^2 



dn (uB,rn) 
1 — TO 
dn 2 (its, to) 



T{uj) = ^5x sin 9 sX /jx 
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FIG. 1: Representative density of states spectrum in the nor- 
mal wire as induced by the proximity effect between two s ++ 
two-band superconductors. Inset shows parameters, and bulk 
gaps were normalized to the Thouless energy Aa = A\/etii- 




FIG. 2: Representative density of states spectrum in the nor- 
mal wire as induced by the proximity effect between two s" 1 
two-band superconductors. The anti-symmetric Fano-like fea- 
tures near the gap edges Aa = Ax/eth may provide a definite 
fingerprint of s +_ -superconductivity. 



Equation ( 10 ) defines m and thus 6q as a function of 
energy oj and together with Eqs. (|3| and ^ it provides a 
complete analytical solution for the Green's function in 
the wire. 

In particular, we can find density of states in the mid- 
dle of the wire as N(e)/Nq = Re[2m(ie) — 1], which 
exhibits very rich structure. Indeed, Fig. [I] shows rep- 
resentative profiles of N(e) in s ++ |n|s ++ junction for 
different choice of parameters. One finds a proximity- 
induced energy gap e g in the spectrum of a wire, which 
scales with the Thouless energy e g ~ Etk- Asymptotic 
analysis near the gap, e — e g -C e g , shows that DOS 
has a square-root singularity N(e) oc ^Jeje a — 1, simi- 
lar to that in a single-band s\n\s junctions.^ N(e) then 
rapidly grows, passes through the maximum and has two 
additional peak-like features at higher energies near the 
superconductive band gaps A>. This picture has to be 
contrasted to the DOS profile in s H \n\s^ junctions 
shown in Fig. [2j The low-energy behavior is similar but 
the energy gap is reduced due to anti-proximity effect in- 
duced by the 7r-shifted band. The conceptual difference 
appears near the band gaps A\ where instead of peaks 
one finds Fano-like antisymmetric features. This impor- 
tant detail is specific for the s" 1 symmetry case and can 
be looked for in the tunneling experiments. 



IV. JOSEPHSON CURRENT 

Phase-sensitive measurements are clearly more chal- 
lenging. Nevertheless, we develop a theory for the 
Josephson effect in mesoscopic s|n|s H circuits with the 
idea that some limits considered here will be useful for 
the future experiments. 

In the presence of a superconducting phase gradient 
in the wire finding an analytical solution of Usadel equa- 
tions represents a difficult technical problem. In the limit 
of the long junctions however L 3> y/D/T calculation 
of the Josephson current simplifies considerably. In this 



case, it is possible to neglect the mutual role of supercon- 
ducting leads and introduce an ansatz for the anomalous 
Green's function F = e^/ 2 sind R + e~ l< ^/ 2 sin8 L , where 
the functions Q R ( L ) satisfy the same sin-Gordon equation 
as in the case of DOS calculations. Solving it separately 
near right (left) boundary for respectively we find 37 



tan[^ L ) (x, w)/4] = B R{L) (w) exp[±(x T l/2)L/& 



(13) 



where we introduced coherence length £ w = y/ D /2uo. 
This approximation conserves the current in the normal 
layer with the exception of the narrow region of the or- 
der £w=t near the boundaries. The two integration coef- 
ficients B r(l) are to be found from the boundary condi- 
tions Eq. (2b) at both interfaces, which can be reduced 



to the algebraic equation: 



AQ a (B a - Bl) - F a {l - QBl + B%) 
= ±2(L/^)(B a +B 3 a ), a^R,,L 



(14) 



where Q and J 7 - functions were defined earlier in Eqs. ( 11 ) 



and (12). With this at hand, we find Josephson current 



in the form 



eI((p)R N = 1287rTsm0^ — B R (uj)B L (uj)e- L/ ^ (15) 



LU>0 



which is applicable in the broad range of temperatures 
£Th <C T < |Aa|. At lowest temperatures T <C em, the 
current-phase relationship in Eq. (15) deviates from be- 



ing simply sinusoidal because a separable approximation 
for ^-function fails to account properly for the proximity- 
induced Thouless gap. Unfortunately, analytical calcula- 
tion of /(</>) is not possible in this limit, however one may 
easily estimate the magnitude of the critical current as 
e/ c i?Ar ~ STh- Furthermore, it is expected that I ((f)) will 
be nonanalytical function near </> — 7r since the proximity 
gap closes at that point while the current is proportional 
to its derivative I(<j)) oc d^Egijp). 
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FIG. 3: (Top) Critical current for s|n|s H junction vs bound- 
ary resistance mismatch r 7 . In the inset, A 3 labels the gap of 
an ordinary superconductor, while Ai,2 are the two gaps of an 
s H superconductor. The other parameters are T = 0.5A S , 
71 = 7s = 5, and L/^t = 2. (Bottom) Critical current 
for s|/|s junction vs boundary resistance mismatch r 7 for 
T = 0.3A S , Ai = 0.5A S , A 2 = 1.5A S , and h = 3A S . 



We plot in Fig. [3] the critical current from Eq. (15) as 
a function of the ratio between the interface barriers for 
each band r 7 = 71/72- It is well known that for the con- 
ventional superconductors, I c decays monotonously with 
r 7 , however, for alnls " 1 junctio n critical current displays 
clear — w switching! 23 * 25 1 28 1 3 ^ This effect is magnified in 
the presence of ferromagnetic layer. Including the ex- 
change field h in Eq. (la I as lo — > u + i/isign(cj), but ig- 



noring spin-flip and spin-orbital scattering, we find from 
the linearized Usadel equations the current in s | | s -1 
junction I(<p) = I c sm(j) with 



eI r R 



N 



(L /^\)(T R T L /g R g L )/cosh(L/^ H ) 
(l + r 2 )tanh(L/£ M )+r^ 

(16) 

where T w = L/2^y/Q R Q L and (j, u = (Qr + Ql)/VGrGl- 
Equation ( |16[ ) is the generalization of the Buzdin for- 
mula^ for the multi-band case. Lower panel of Fig. [3] 
shows enhanced 0-7r oscillations of the critical current as 
a function of r 7 , which displays two zero points. Such a 
peculiar feature is due to the combination of a ferromag- 
net and s~* superconductor. 

In the Josephson junction with extremely low bar- 
rier transparency when 7^ S> 1, one can circumvent the 
need of solving Usadel equation in the wire since cur- 
rent is largely determined by the interface. Supercon- 
ductive phase 4> changes discontinuously at the barri- 
ers and stays nearly zero within the interior of the wire 
while Green's function phase 9 is approximately constant. 
Since J w oc 7 -1 <c 1, then to the leading order, one can 
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FIG. 4: Josephson current-phase relationship for short 
s|?i|s + ~ junction for gaps Ai = 0.5A S , A2 = 1.5A S , and 
interface parameters j a = 2, 71 = 2.4, 72 = 1.8. 



set X b = in Eq. (2a I and d x = in Eq. (2b). These 



approximations allow to find the spectral current, 



J u = T sin ( 



T cos 2 



-1/2 



(17) 



where we assumed s" 1 \n\s~ 1 geometry. For the equal 
gaps case Ai = A2 = A s , above the spectral current J u 
leads to the Josephson current-phase relationship from 
Eq. @: 



eI{(j>)R N 



M 2 (l + r 7 ) 



27i 



A, sin < 



tanh(e/2T)(fe 



fj,A s cos * ^A 2 - e 2 J e 1 - A 2 cos 



(18) 



where r 7 = 71/72 and \x = (1 — r 7 )/(l+r 7 ). Interestingly, 
in the zero-temperature limit, even the remaining energy 
integral can be completed in the closed form, such that 
we find a Josephson current 



eI(4>)R N 



(72 " 7i) 2 A s sin0 
7172(71 +72) 



K 



72 - 7i 
72 + 7i 



cos 



2 
(19) 

where K(x) is the complete elliptic integral of the first 
kind. In the completely symmetric case (with equal 
transparencies for both bands), the Josephson current 
vanishes, since the 7r-shifted bands drive it in the oppo- 
site directions. However, including interband scattering 
effects may result in additional nonvanishing contribu- 
tions P3 

Analytical results are also possible for arbitrary trans- 
parencies but for the short junctions when L <C £,ui=A x ■ 
In this case, Usadel equations are dominated by the gra- 
dient terms. Despite the fact that they remain coupled 
and nonlinear, all integrations can be completed in the 
closed formPS By using the first integral of Eq. ( lb ) , 
J u = sin 2 9d x x and excluding d x x from Eq. (la), one 
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finds 



where 



2sin 4 6> 



(20) 



This nonlinear differential equation is solved by 

cos[6(x,w)] — cos9oCos[J LJ (x ~ x )/ sinOo]. (21) 

Knowing 9{x,oS), one can now calculate the second inte- 
gral of Eq. ([Ibl , 



X - Xo = Ju. 



x dx 



sin 2 9 ' 



(22) 



which reads 

sin9 tan[x(x,uj) - xo] = ta,n[J u (x - x )/sm9 ]. (23) 

Having found explicit solutions for the Green's functions, 
the boundary problem for the integration coefficients can 
be reduced to solving three algebraic equations: 

J u = J~ l sin 9 B sm( XB + 0/2) , (24) 

J u = J^r sin 6 B sin(0/2 - xb), (25) 
Tl cos 9 b cos(xb + 0/2) - Q L sin 0b 

= F r cos9bCOs(xb - 0/2) - GRsm9 B - (26) 

These expressions finally lead us to the Josephson 
current-phase relationship in the form 



eI{<t>)R N = 8ttT 



A(4>) sin 



Ft 



A 2 W 



(Or - Ql? 
(Tr - Tl? 



(27) 



-4(0) 



{Fr-Fl? . 2 

(T R + F L y sm 2 



-1/2 



(28) 



A representative feature of Eq. (27) is that 7(0) switches 

7r as shown in Fig. [4] 
junction has 



its sign in between = and 
This implies that the free energy of s|n|s + " 
two minima and such junction may be used as the phase 
inverted in superconducting digital circuits. Such fea- 
ture, however, is not unique for s H superconductors and 
can be realized in other complex hybrid circuits with or- 
dinary materials.^ 

In summary, we have studied the density of states and 
Josephson current in mesoscopic circuits with unconven- 
tional s H superconductors. We find that tunneling spec- 
tra have distinct fingerprints of the sign changing sym- 
metry of the underlying superconductive order parameter 
induced by the proximity effect. Furthermore, the critical 
current exhibits a robust n junction even in the absence 
of the ferromagnetic layer. The Josephson current-phase 
relationship itself is not indicative of symmetry due 
to the sensitivity to parameters defining the junction. 
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ing on the paper. This work was supported by Michigan 
State University. 
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